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The R-matrix method has been used for theoretical calculation of electron collision 
with atoms and molecules for long years. The method was also formulated to treat 
photoionization process, however, its application has been mostly limited to photoion- 
ization of atoms. In this work, we implement the R-matrix method to treat molecular 
photoionization problem based on the UK R-matrix codes. This method can be used 
for diatomic as well as polyatomic molecules, with multi-configurational description 
for electronic states of both target neutral molecule and product molecular ion. Test 
calculations were performed for valence electron photoionization of nitrogen (N2) as 
well as nitric oxide (NO) molecules. Calculated photoionization cross sections and 
asymmetry parameters agree reasonably well with the available experimental results, 
suggesting usefulness of the method for molecular photoionization. 
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I. INTRODUCTION 



Photoionization of molecule has been studied for long year because of its importance in 
modeling and understanding upper atmosphere, interstellar clouds and industrial plasma. 
Photoionization process has also been used to understand nature of molecular electronic 
states as well as dynamics of nuclei on excited electronic stated In addition to the con- 
ventional measurement of integral cross section and photoelectron asymmetry parameter 
for randomly oriented molecule, recent developments of experimental technique have made 
it possible to extract photoelectron angular distribution from oriented molecule^, perform 
accurate measurement of inner-shell photoionization cross section 3 . 

Several theoretical methods have been developed for treating photoionization of molecule, 
based on the Schwinger variational principle^, random phase approximation^, time- 
dependent density functional theory^, complex basis function 8,9 , Stieltjes imaging method^, 
and so on. These methods have successfully reproduced or predicted experimental results 
for photoionization of small to medium sized molecules in the gas phase. 

The R-matrix method, originally developed for nuclear reaction^, has been applied for 
accurate calculation of cross section for electron collision with atoms, molecules and ions 
with great success^— . As the other electron scattering theories such as the Schwinger 
variational method^ and the Kohn variational method^, the R-matrix method was also 
formulated to treat photoionization of atoms and molecules-^. In addition to the single- 
photon ionization process, the method was also modified to treat multi-photon ionization 
process^ 9 - 1 ^. There have been many applications of the R-matrix method to photoionization 
of atoms^^ 8 -, however, its application to molecular photoionization has been fairly limited, 
with only a few reports on hydrogen molecule^ 1 ^. The ab initio R-matrix method can 
reproduce the existing experimental cross sections very well, as has been demonstrated in 
the previous electron-molecule scattering calculations^ 3 - 1 ^. Thus, the method is expected 
to work efficiently in molecular photoionization problem as well. Introduction of a differ- 
ent method, the R-matrix method, may be useful when various theoretical methods are 
compared with experiment. In the R-matrix treatment of photoionization, it is possible to 
represent electronic states of target molecular ion by multi-configurational wavef unctions, 
as it has been in the R-matrix calculations on electron molecule collisions^ - — This 
multi-configurational representation of wavefunction may be good for description of satellite 
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or resonance states with multi-configurational character, as well as introducing correlation 
in the initial and final states. 

In this work, we implement the R-matrix method of molecular photoionization based 
on the work of Burke and Taylor 18 originally proposed for atomic photoionization. Our 
procedure is roughly divided into two stages. In the first stage, we perform R-matrix cal- 
culations for electron collision with molecular ion, and obtain eigenvalues and eigenvectors 
of the R-matrix eigenstates by diagonalizing the electronic Hamiltonian inside the R-matrix 
sphere. In the second stage, the initial and final state wavefunctions of photoionization are 
constructed as linear combinations of the R-matrix eigenstates obtained in the first stage, 
then the transition dipole moments between these initial and final state are obtained for cal- 
culation of photoionization cross section and asymmetry parameter. For calculation in the 
first stage, we use the polyatomic version of the UK R-matrix codes developed by Morgan 
et al.— The UK R-matrix codes has been successfully applied to many electron-molecule 
scattering problems in the past. In this work, some modifications are made to the UK 
R-matrix codes, e.g., to enable transition dipole moment calculation between the R-matrix 
eigenstates. In order to verify reliability and accuracy of the R-matrix method for molecular 
photoionization, cross sections and asymmetry parameters are calculated for valence elec- 
tron photoionization of nitrogen (N2) and nitric oxide (NO) molecules randomly oriented 
in the gas phase, and the results are compared with existing experimental and theoretical 
data. N2 has been benchmark molecule for various theories of photoionization, and will be 
also good for the first test calculation of the R-matrix method for molecular photoionization. 
In contrast to N2 molecule, number of calculations on photoionization of NO molecule is 
limited^. We selected NO molecule as the second test case, to see capability of the method 
to treat this open-shell molecule with II ground state. 

II. THEORETICAL METHOD 

In this section, we describe the procedure to obtain molecular photoionization cross sec- 
tion and asymmetry parameter using the R-matrix method. Our discussion mainly follows 
Chandra^^. and Burke and Taylor 18 . Many equations in this section already appear in their 
papers, however, they are shown here for reader's convenience. 

Differential cross section for photoionization of molecule randomly oriented in space can 
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be written as,— ^ 
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where fc' represents direction of a photoelectron with wavenumber in laboratory frame, 
E r is photon energy, is a Legendre polynomial, 9' is the angle between the direction of 
photoelectron and the polarization vector of the incident photon. Here we assume that the 
photon beam is linearly polarized. The expansion coefficient Ai(k) is represented as 
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Here e\ r is the unit vector describing polarization A r of the incident photon in molecular 
frame, n e is number of electrons in the system with r s being coordinates of the sth electron. 
In this work, we employ dipole length approximation. Thus, the proportional coefficient A 
becomes^ (47ra 3 £'^/3e 4 ) 1 ^ 2 with a being the fine structure constant. As is well known, the 
differential photoionization cross section in Eq. ([T]) can be simplified in the form, 
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with the integrated cross section, 
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Following Burke and Taylor—, we evaluate the transition dipole matrix elements in Eq. ([3]) 
using the R-matrix method^— for electron-molecule collision. In this method, configuration 



space is divided by two regions according to the distance r„ e of the scattering electron, i.e., 
photoelectron, and the center of mass of the molecular ion having (n e — 1) electrons. In 
the inner region, defined by condition r Ue < a, the n e -electron wave functions of the total 
system are represented by (n e — l)-electron wave functions of molecular ion augmented by 
diffuse functions which describe a scattering electron, 

+ ^2x q (x 1 ...x ne )bq, (7) 

where A is an antisymmetrization operator, Xj represents spacial and spin coordinates of the 
ith electron, 0, are the eigenstate of the total spin and its z component constructed from 
the (n e — l)-electron wave functions of the molecular ion 0, and the spin function of the 
scattering electron, Uj are continuum orbitals representing wave functions of the scattering 
electron, and X q are bound n e electron wave functions composed of the target molecular 
orbitals and the extra target virtual orbitals, the coefficients ajj and b q are determined by 
diagonalization of H ne + L ne where H ne is the electronic Hamiltonian and L ne is a Bloch 
operator— 1 ^ accounting for surface term. In the outer region r ne > a, the problem is 
reduced to single electron scattering, ignoring exchange of the scattering electron with the 
electrons of the molecular ion. Interaction of the scattering electron and the molecular ion 
is considered through static multipolar interaction potentials which introduce inter-channel 
couplings. The R-matrix eigenstates obtained by diagonalization of the Hamiltonian in the 
inner region are converted to the R-matrix Rij at the boundary r ne = a as, 
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where is eigenvalue of E is the energy of the total system, and the boundary ampli- 
tudes Wik is given by 

w ik (r) = ($iY limi \^ k ). (9) 

For photoionization problem, the final states ^{~ mf i n Eq. ([3]) are expanded by the R-matrix 
eigenstates $£: as^, 
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Here ^{~ mf is denoted as for simplicity. Using the relation^^ which holds for the 
R-matrix amplitudes and the final state wave functions, the expansion coefficients can be 



written as 



where Ef is the energy of the final state, which is equal to E^ + E r with the initial state 
energy E{ and the photon energy E r . The parameter b is related to the logarithmic derivative 
of the scattering electron wavefunction at the R-matrix boundary and is set to zero in this 
work. The radial function yjf is defined as 
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We have to determine the radial functions and their derivatives in Eq. (TTTT) to evaluate the 
final state wave functions. In the outer region, the radial functions yjf satisfy the differential 
equations^^, 
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where n is the number of channels considered in the R-matrix model, z is the net charge 
of the molecular ion, kj is the wavenumber of the electron in the jth channel and Vjk (r) 
represents multipole potential.^ When n a channels are open, we have n + n a independent 
solutions of Eq. (f!3l) with the asymptotic boundary conditions, 
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with rji = —z/ki. The radial functions yjf are expanded by these n+n a independent solutions 
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By inserting Eq. ( fl6l) into the R-matrix relation— 1 ^, 
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and using the ingoing wave asymptotic conditions of yf, 
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with K and being K-matrix and its elements, we obtain (n + n a ) linear equations for 
(n + n a ) unknown coefficients x k . The expansion coefficients x k can be determined by solving 
these equations, and as a result, we can evaluate the final state wavefunctions <3>J using Eqs. 
(HUD and (HU. 

In this work, the initial state wave function is also expanded by the R-matrix eigen- 
states as described in Burke and Taylor—, 

^i = J2 A ki^k, (19) 

k 

where the coefficient A ki is given by Eq. ffTTj) with Ef and yjf replaced by the energy Ei 
and the radial function of the initial state. The radial functions yji are expanded by the 
independent solutions of the differential equations (IT5|) as, 

n 

y j i = ^2 / v jk x k . (20) 
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Since we are treating bound initial state, all channels are closed, n a = 0, and all independent 
solutions decay to zero as r approaches infinity. By substituting Eq. (|2"U|) into the R-matrix 
relation of Eq. ( fl7|) . n equations for n unknown coefficients x k are obtained. These equations 
are only solved at discrete energies of Ei, where the lowest of them corresponds to the ground 
state energy of the neutral molecule. 

III. APPLICATION OF THE METHOD TO PHOTOIONIZATION OF N 2 
AND NO MOLECULES 

A. Detail of the calculations 

Cross sections and asymmetry parameters for valence electron photoionization of N 2 and 
NO molecules were calculated based on the method described in the previous section. The 
R-matrix eigenstates and amplitudes were obtained by electron - and NO + scattering 
calculations using a modified version of the polyatomic programs in the UK molecular R- 
matrix codes^. We used the fixed nuclei approximation with internuclear distances 2.068 ao 
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for N2 and 2.175 ao for NO, which are the equilibrium values of the ground electronic states 
of N 2 and NO. 

For the electron - Nj R-matrix scattering calculation, we employed two R-matrix models; 
a single-channel SCF target model and a multi-channel CI target model including 80 N^ 
electronic states. The cc-pVTZ atomic basis set^i was used to describe molecular orbitals 
in both models. In the CI target model, full valence complete active space was employed to 
generate configuration state functions for the N^~ electronic states. The molecular orbitals 
(MOs) in this active space were obtained by the state-averaged complete active space self 
consistent field (CASSCF) 32 ' 33 calculation using molpro program package 3 ^. The state- 
averaging was performed over the lower 28 electronic states of N^~ ion, then the CASCI 
wave functions of the remaining 52 electronic states were constructed from this CASSCF 
molecular orbital set. The CASSCF ionization potentials of the lower 9 N^ electronic states 
are listed in Table [U where the N2 ground state energy was obtained by the electron - 
N^~ scattering calculation, as described at the end of the last section. In addition to these 
CASSCF MOs for the N^ valence electronic states, we included 3 extra virtual orbitals for 
each irreducible representation of the D 2 h symmetry. In order to represent the scattering 
electron, we included diffuse Gaussian functions up to I = 5, with 12 functions for / = and 
1, 8 functions for / = 2 and 3, and 5 functions for I = 4 and 5. Exponents of these diffuse 
Gaussians were obtained by the GTOBAS program 3 ^ in the UK R-matrix codes. To prepare 
an orthogonal MO set used in Eq. (I7j), these Gaussian functions were orthogonalized against 
the valence and extra virtual MOs obtained by the CASSCF calculation. The procedure for 
construction of the 14-electron configurations in Eq. ([7]) is almost the same as we did in the 
previous work of electron - N 2 scattering calculation^, except the number of electron in the 
system, and is not repeated here. Radius of the R-matrix sphere was chosen to be 10 ao 
in our calculations. The R-matrix calculations were performed for the singlet A g , Bi u , B2 U 
and B^ u symmetries of the e+N^ system, where the A g result was used for the ground state 
of the neutral molecule whereas the other symmetries were used to construct the final state 
wavefunctions. 

Detail of the electron - NO + R-matrix scattering calculation is similar to the electron 
- N^ scatterings. In this case, we performed a single-channel SCF target calculation and 
a multi-channel CI target calculation with 60 NO + electronic states. The state-averaged 
CASSCF calculation with cc-pVTZ basis set was performed for the lowest 16 electronic 
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states of NO + , then the CASCI wave functions were constructed for the other states. The 
ionization potentials of the lowest 8 NO + electronic states are shown in Table [Til Radius of 
the R-matrix sphere and diffuse Gaussian functions for the scattering electron are the same 
as in the e-N^ calculation. In this e - NO + calculation, we included 4 extra virtual orbitals 
for each irreducible representation of the C2 V symmetry. The R-matrix calculations were 
performed for the doublet Ai, A<i and B\ symmetries of the e+NO + system, where the B\ 
result was used to describe both the ground state of the neutral molecule and the final state 
wavef unctions, whereas the other symmetries were only used for the final states. 

The wave functions of the initial N 2 X 1 E+ and NOX 2 II states were expanded by the 
R-matrix eigenstates taken from the e - N^ and e - NO + R-matrix scattering calculations 
as described at the end of the last section, and their expansion coefficients were obtained 
by the BOUND module of the UK R-matrix codes. It was found that the wave functions 
of these N2 and NO states are dominated by the lowest eigenvalue R-matrix eigenstates, 
i.e., expansion coefficients of the lowest eigenvalue states are almost unity whereas they 
are less than 10 -4 for the other R-matrix eigenstates. Thus, we just substituted the initial 
state wave functions in Eq. (j3J) by these lowest eigenvalue R-matrix eigenstates. Note that 
the above procedure means that the ground state wavefunction of the neutral molecule is 
described by the MOs of the molecular ion. This is approximate method for description of the 
ground state wavef unctions, however, it greatly simplifies evaluation of the transition dipole 
moments, since the same set of MOs is used in the initial and the final state wavefunctions. 

In order to obtain the transition dipole matrix elements of Eq. ([3]), we evaluated the 
dipole matrix elements between the R-matrix eigenstates. Then the expansion coefficients 
for the final state wave functions were calculated using Eqs. (fill and ( fl6l) . In principle, 
accurate solutions of differential equations (ITS]) with proper boundary conditions ( fT4T) have 
to be used for Eq. (Tl6|) . Such solutions are usually prepared by the asymptotic expansion^&^Z 
and inward integration to the matching point r = a. However, the asymptotic expansion 
is not so accurate for small r and near the ionization thresholds, where we have 80 such 
thresholds for N 2 photoionization and 60 for NO molecule in our R-matrix models. Also, 
inward integration is unstable due to exponential growing of closed channel components. 
In this work, we approximated Vij by the Coulomb functions ignoring multipole potentials 
of the molecular ions. Based on the calculated transition dipole moments, photoionization 
cross sections and asymmetry parameters were evaluated using Eqs. (j2J), (jSJ) and flH]). 
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B. Results and Discussion 



1. Ionization Potentials 

In Table HJ the CASSCF ionization energies of N 2 molecule are shown for the lower 
ion states with the experimental values^. Our results agree reasonably well with the exper- 
imental ionization energies, though ours are 0.5 - 1.0 eV higher. The ionization potentials of 
NO molecule are also shown for the lowest 8 NO + ion states in Table [Til Since the equilib- 
rium bond distances of some excited states are longer than that of the neutral NO ground 
stated, energetic order of the A 1 !!, A fl T,~ and W l /S. states in our results is different from 
that in the adiabatic experimental results. Nevertheless, our IPs compare reasonably well 
with the experiment, with deviations less than 1 eV. Note that the IPs in Tables HI and HT1 
were extracted from the ionic and neutral energies used in the R-matrix photoionization 
calculations. Since we did not attempt to adjust these energies to the experimental IPs, the 
ionization thresholds in our cross sections and asymmetry parameters are shifted relative to 
the correct experimental values. 

2. N 2 

In Fig. [TJ cross section and asymmetry parameter for photoionization of the N 2 X 1 E+ 
state leading to the N 2 f X 2 E+ state are shown. In the figure, our CASSCF results and 
SCF results are compared with the previous theoretical results of Montuoro and Moccia 6 
obtained by the K-matrix method with interacting channels random phase approximation 
and Stratmann et al.— obtained by the multi-channel Schwinger variational method. Also, 
the experimental cross sections of Hamnett et al.— , Samson et al.— , Plummer et al 4 ^., and 
the asymmetry parameters of Marr et al. 44 and Southworth et al.— are included in the figure. 
Our SCF cross section is smooth and has a large broad peak around 30 eV, originated from 
the a* shape resonance. In contrast, our CASSCF multi-channel cross section has numerous 
sharp peaks originated from two- or many-electron excited resonances as well as Rydberg 
resonances associated with the excited electronic states of N^ ion. Overall shape of the cross 
section is roughly similar to the SCF result. The CASSCF cross section rises at low energy 
region around 15-20 eV, which is consistent with the experimental result of Hamnett et al. 
and Samson et al. and the previous theoretical cross section of Montuoro and Moccia. Below 
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30 eV, our CASSCF cross section agrees well, except for the presence of the sharp peaks, with 
the experimental results. However, it overestimates the experimental values above 30 eV. 
The SCF cross section agrees rather well with the experimental results in this high energy 
region. Similar to the calculated cross sections, the CASSCF asymmetry parameter /3 has 
numerous sharp peaks while the SCF asymmetry parameter is very smooth. Compared to the 
SCF result, the CASSCF result agrees better with the experimental asymmetry parameters 
in all energy range in the figure. At low energy region below 24 eV, the CASSCF R-matrix 
calculation overestimates the experimental asymmetry parameters, partly due to existence 
of resonances. The cross section and asymmetry parameter of Stratmann et al. are available 
over limited energy range of 19 - 26 eV. In this energy region, our CASSCF results are very 
similar to their results. 

Figure |2] shows photoionization cross section and asymmetry parameter for the N^~v4 2 n 
state, where our CASSCF results are compared with the previous theoretical results of 
Stratmann et al.— and the experimental cross sections of Hamnett et al.— , Samson et al.— , 
Plummer et al.— , and the experimental asymmetry parameter of Marr et al.— Shape of our 
CASSCF cross section is roughly similar to the experimental results. However, magnitude of 
our cross section is slightly larger than the experimental cross sections above 25 eV. In case 
of the asymmetry parameter, our result agrees very well, except for the presence of the small 
peaks, with the experimental result of Marr et al. The previous multi-channel Schwinger 
results of Stratmann et al. are available between 19 and 26 eV. In this energy region, the 
shape of our cross section is very similar to their result, although magnitude of cross section 
is larger than theirs. For the asymmetry parameter, our result almost coincides with the 
result of Stratmann et al. 

In Fig. |3l cross section and asymmetry parameter for photoionization leading to the 
N^~I? 2 £+ state are shown with the available experimental results and the previous theo- 
retical results of Stratmann et al.— As in the photoionization to the N^~A 2 n state, our 
cross section is slightly larger than the result of Stratmann et al., however, the asymmetry 
parameter is very similar to each other. Compared to the experimental results, our cross 
section looks larger near the threshold, yet, agreement is better at higher energies. On aver- 
age, our asymmetry parameter is roughly similar to the experimental results, however, the 
experimental asymmetry parameter has a large dip around 30 eV which does not exist in 
our results. In the previous studies&^i, a coupling between the A 2 £+ and B 2 T>+ channels 
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has been suggested to cause this large dip. Although such a coupling is included in our 
model, our asymmetry parameter is rather flat between 20 and 45 eV, with slight decrease 
of magnitude around 35 eV. 

In general, partial photoionization cross sections for the higher ion states are very 
small compared to those for the N^X 2 E+, A 2 U and B 2 T^ states. However, the photoioniza- 
tion cross section leading to the 2 2 E+ state is relatively large. The calculated cross section 
and asymmetry parameter are shown in Fig. HI The ionization energy of this N^2 2 E+ state 
is about 29.65 eV, which is close to the threshold of the "Z" state mentioned in Hamnett 
et al.— Their cross section for the "Z" state is also shown in Fig. HJ Although agreement 
is not so good, the magnitudes of the cross sections are roughly similar to each other. The 
calculated asymmetry parameter drops from 1.5 at the threshold to -0.5 near 35 eV, then 
it increases to 0.5 at 45 eV. This behaviour differs from the asymmetry parameters for the 
N^"A, A and B states, probably reflecting the difference in the main electronic configurations 
of the ionic states. 

3. NO 

In Fig. [51 cross section and asymmetry parameter for photoionization of the NOX 2 IT 
leading to the NO+X 1 ^" 1 " state are shown with the experimental results of Southworth 
et al. 46 and Iida et al.— , and the previous theoretical results of Stratmann et al.— Our 
CASSCF cross section has numerous sharp peaks between 13 and 18 eV, mostly originated 
from Rydberg resonances associated with the excited electronic states of NO + . In the other 
energy region, the magnitude of the cross section is nearly constant value of about 5 Mb. 
Compared to the CASSCF result, the SCF cross section is very smooth and flat. Although 
a slight increase of the SCF cross section exists around 16 eV, the height of the peak is 
much smaller than the CASSCF cross section around this energy region. The shape of our 
CASSCF cross section is very similar to the result of Stratmann et al, especially above 14 
eV. Between 12.5 and 14eV, our cross section has several sharp peaks as in the higher energy 
region. In contrast, there is no such peak in the cross section of Stratmann et al. below 
14 eV. Below 25 eV, our CASSCF cross section agrees well, except for the presence of the 
sharp peaks, with the experimental cross section of Southworth et al. Our result, however, 
slightly underestimates experimental cross section above 25 eV. Compared to the result of 
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Iida et al, the magnitude of our CASSCF cross section is generally smaller except in the 
high energy region around 35-40 eV. Enhancement of cross section is observed below 20 eV 
in the experimental results of Southworth et al. and Iida et al., which is roughly reproduced 
by our CASSCF R-matrix calculation. However, this feature is not well captured in the SCF 
R-matrix model. The shape of our CASSCF asymmetry parameter is roughly similar to the 
result of Stratmann et al. The asymmetry parameter of Stratmann et al. is rather steeper 
above 18eV, and agrees better with the experimental result of Southworth et al. Our SCF 
asymmetry parameter has a much smoother profile compared to the CASSCF result. Its tilt 
is close to the results of Stratmann et al. and Southworth et al., though the magnitude is 
generally larger. 

Figure [6] shows cross section and asymmetry parameter for ionization to the NO + & 3 IT 
state. The experimental results of Southworth et al.— are also included in the figure. Near 
the ionization threshold, our cross section has many peaks as in the cross section of the 
NO + X 1 S + state. The magnitude of the cross section increases from 5 Mb to 7 Mb as 
photon energy increases from the threshold to 27 eV. Then it decreases to about 3 Mb at 
40 eV. Although overall shape of the cross section is similar to the experimental result, the 
broad peak in our result is located 3-4 eV higher in energy compared to the peak in the 
experimental cross section. Our asymmetry parameter decreases from 1.5 to 0.4 as energy 
increases from the threshold to 40 eV. This behaviour is roughly similar to the experimental 
result, however, our asymmetry parameter has a dip around 20-25 eV whereas there is a 
small bump in this energy region in the experimental result. 

In Fig. [71 cross section and asymmetry parameter for photoionization to the NO + A 1 II 
state are shown with the previous theoretical results of Stratmann et al.— Our cross section is 
generally larger than their result above 22 eV. The shapes of the cross sections resemble each 
other, having a broad peak around 25 eV. The asymmetry parameter gradually decreases 
from 1.5 at the threshold to 0.15 at 40 eV. Our result has a dip around 23 eV which does 
not exist in the asymmetry parameter of Stratmann et al., however, the other details are 
similar. 

In Fig. [El photoionization cross sections and asymmetry parameters for the NO + w 3 A, 
a 3 S + , A' 1 ^ , 5' 3 £~ and W /1 A states are shown with the averaged experimental asymmetry 
parameter of Southworth et al.— The cross sections for the a, A', b' and W states are 
relatively flat above the thresholds and are similar to each other. The cross section for the 

13 



w state decreases slowly from the threshold to 40 eV, with a drop around 20-23 eV. The 
magnitude of cross section is the largest for the w 3 A state with average value of about 5 Mb. 
The magnitude of the cross sections for the a 3 £ + , 6 /3 S~ and W 1 A states is about 2.0-2.5 
Mb, and the A fl Tj~ state has the smallest magnitude of about 1.0 Mb. The asymmetry 
parameters of these states have remarkably similar shapes and magnitudes, starting from 
-0.8 at the thresholds and increasing monotonically to 1.2 around 40 eV. Our results agree 
well with the experimental averaged asymmetry parameter of Southworth et al. 

4- Discussion 

The photoionization cross sections and asymmetry parameters calculated in this work 
have numerous sharp narrow peaks. They are originated from the Rydberg resonances 
associated with the excited electronic states of the molecular ions as well as two- or many- 
electron excited resonances. This can be checked by changing number of excited ionic states 
included in the R-matrix calculation or by modification of the X q terms in Eq. (TH). These 
sharp resonances are characteristic of multi-channel method such as present CASSCF multi- 
channel R-matrix calculation and the Schwinger multi-channel method. Since they are not 
described in the single-channel SCF target R-matrix calculation, the SCF cross section 
and asymmetry parameter in Figs [1] and [5] are very smooth. These numerous peaks are 
rarely seen in the previous theoretical calculations because many of them used single-channel 
method. The cross sections and asymmetry parameters of Stratmann et al.— , calculated by 
the Schwinger multi-channel method, are very similar to our CASSCF results, including 
location of several peaks. For NO photoionization, some energy regions exist where our 
results have many narrow peaks whereas the results of Stratmann et al. are smooth. This 
discrepancy is attributed to difference of NO + electronic states included in the calculations. 
For example, Stratmann et al. did not include NO + ionic states with (tt)^ 1 configuration, 
however, we put these electronic states in the R-matrix model. These narrow peaks in our 
results are not observed in the previous experimental results. In this work, we employed the 
fixed nuclei approximation. However, if we include the effect of vibrational motion by using 
the adiabatic averaging method or the non-adiabatic R-matrix method, these narrow sharp 
peaks may be averaged and look less prominent. 

As can be seen from the figures, the shapes of the photoionization asymmetry parameters 
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are closely related to the main electronic configuration of the molecular ion, in other words, 
the molecular orbital from which ionization occurs. When the molecular ion has a (cr) _1 con- 
figuration, e.g., N^~X 2 £+ and NO + 6 3 IT states, the asymmetry parameter is generally 0.5-1.0 
near the threshold and remains nearly constant or decreases slightly with energy. In contrast, 
when the configuration is (7r) _1 type, e.g., N^t4 2 II u and NO + X x S + states, the asymmetry 
parameter is negative or zero near the threshold and increases rapidly with photon energy. 
These behaviour have been discussed by several authors including Southworth et al.— and 
Thiel^, and are interpreted in terms of angular momentum components of photoelectrons 
near threshold. 

In our model, angular momentum of photoelectron was considered up to I = 5. Photoion- 
ization cross sections and asymmetry parameters obtained by maximum angular momentum 
of I = 3 and 5 are very similar to each other, which suggests that we do not have to increase 
I more than 5 in this study. Concerning number of excited electronic states of molecular ion, 
we included 80 and 60 electronic states of N^" and NO + ions, respectively. This inclusion of 
many excited ionic states is necessary to obtain converged cross sections and asymmetry pa- 
rameters as well as to suppress pseudo resonances. The numbers of states in our models are 
much larger than those in the previous Schwinger multi-channel calculations of Stratmann 
et al.— ^ In the asymptotic region, we considered only the Coulomb potential in this work. 
In order to inspect the effect of dipole and quadrupole potentials on resonance position, 
we calculated elastic cross sections of electron - and NO + collisions with and without 
the multipole potentials, using CASSCF multi-channel R-matirx models. As in the pho- 
toionization cross sections shown in the figures, the elastic cross sections of the electron-ion 
collisions contain many sharp resonance peaks. As far as we have checked, the effect of the 
multipole potentials on resonance position in these elastic cross sections is not significant, 
less than 0.01 eV. This observation suggests that arrangement of resonances in photoioniza- 
tion cross section does not change much by introduction of the multipole potentials, at least 
for N 2 and NO molecules. However, when accurate assignment of Rydberg resonances and 
determination of quantum defects are of interest, consideration of dipole and quadrupole 
potentials will be important. 

In this work, we used the SA-CASSCF MOs of the molecular ion, or NO + , to describe 
the ground electronic state of the neutral molecule, N 2 AT 2 £+ or NOX 2 II. Using a single set 
of MOs for both the initial and final states makes evaluation of transition dipole matrix 
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elements much easier and faster than using non-orthogonal MO sets optimized for neutral 
and ionic electronic states separately. Our SA-CASSCF MOs are not optimized for neutral 
molecule, thus, the full-valence CASCI energy of the neutral molecule obtained by these 
ionic MOs is generally higher than the energy obtained by neutral MOs. Since we have 
extracted the energy of the neutral molecule by the R-matrix scattering calculation, the 
ground state wavefunction of the neutral molecule is constructed from the configuration 
with the valence MOs as well as the additional configurations which contain diffuse orbitals, 
as described by Burke and Taylor— and explained at the end of the last section. Inclusion 
of these additional configurations generally lowers the ground state energy compared to the 
full-valence CASCI energy. For valence electron photoionization, the effects of these ionic 
MOs and the additional configurations appear to compensate each other, since deviations of 
the calculated IPs from the experimental values are not significant as shown in Table [J and 
HT1 However, in case of inner-shell ionization, this compensation does not work well, e.g., 
the ionization threshold for the N^"(la 9 ) _1 state is about 4 eV lower than the experimental 
value in our preliminary calculation. Thus, we may need to use non-orthogonal MO sets 
to describe neutral and ionic states, when the R-matrix method is applied to inner-shell 
photoionization. 

Comparison of our calculation with the previous experimental and theoretical works 
demonstrates that the R-matrix method is well suited to treat valence electron photoion- 
ization of closed-shell as well as open-shell molecules such as N 2 and NO. We implemented 
the method as simple as possible. So, effects of the velocity form dipole moment, multi- 
pole potentials, vibrational motion etc. are not discussed or considered in this paper. In 
the future, we will study these issues along with application of the method to inner-shell 
molecular photoionization. 

IV. SUMMARY 

In this work, we implemented the R-matrix method to treat molecular photoionization 
based on the procedure described by Burke and Taylor—. For the inner region calculation, 
the polyatomic version of the UK R-matrix codes was used with some modifications. Final 
state wavefunctions of molecular photoionization were represented by linear combinations 
of the R-matrix eigenstates obtained by diagonalization of the electronic Hamiltonian inside 
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the inner region. Then, transition dipole matrix elements between the initial and the final 
states were calculated for evaluation of observables such as photoionization cross section 
and asymmetry parameter. As test calculations, we applied this method to valence electron 
photoionization of N 2 and NO molecules. In our R-matrix models, the electronic states of 
the molecular ions, and NO + , were represented by full- valence complete active space CI 
wavefunctions with cc-pVTZ basis set. Calculated photoionization cross sections and asym- 
metry parameters have many narrow peaks originated from two- or many-electron excited 
resonances as well as Rydberg resonances associated with the excited electronic states of N^" 
or NO + , which are very similar to the previous theoretical results of Stratmann et al.— ^ 
obtained by the multi-channel Schwinger calculations. Overall shapes of the calculated cross 
sections and asymmetry parameters agree reasonably well, except for the presence of the 
sharp peaks, with the available experimental results. 
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LIST OF FIGURES 



1 Photoionization cross section (a) and asymmetry parameter (b) for ionization 
of the N^~X 2 £+ state. The cross sections are shown in unit of mega barn 
(Mb), equal to 10~ 22 m 2 . Our CASSCF target results and SCF target results 
are represented as thick full lines and thick dashed lines, respectively. The 
previous theoretical results of Montuoro and Moccia£ and Stratmann et al.— 
are shown as thin dotted and thin dashed lines, respectively. Experimental 
data included in the panel (a) are taken from Hamnett et al.— , Samson et 
al.— and Plummer et al^. Experimental data in the panel (b) are taken from 
Marr et al.— and Southworth et al.— |22 

2 Photoionization cross section (a) and asymmetry parameter (b) for ionization 

of the N^t4 2 I1 state. The other details are the same as in Fig. [TJ [23 

3 Photoionization cross section (a) and asymmetry parameter (b) for ionization 

of the N^~.B 2 £+ state. The other details are the same as in Fig. [TJ |24 

4 Photoionization cross section (a) and asymmetry parameter (b) for ionization 
of the N^~2 2 E+ state. The symbol represents the cross section for the "Z" state 

of in Hamnett et al.— The other details are the same as in Fig. [TJ |25 

5 Photoionization cross section (a) and asymmetry parameter (b) for ionization 
of the NO + X 1 S + state. Our CASSCF target results and SCF target results 
are represented as thick full lines and thick dashed lines, respectively. The 
previous theoretical results of Stratmann et al.— are shown as thin dashed 
lines. Experimental data in the figure are taken from Southworth et al.— and 

Iida et al& I26 

6 Photoionization cross section (a) and asymmetry parameter (b) for ionization 

of the NO + 6 3 n state. The other details are the same as in Fig. [5j |27] 

7 Photoionization cross section (a) and asymmetry parameter (b) for ionization 

of the NO+A 1 !! state. The other details are the same as in Fig. [5j [28| 
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Photoionization cross section (a) and asymmetry parameter (b) for ionization 
of the NO+w 3 A, a 3 S+, A n T,~ , 6 ,3 £- and W X A states, obtained by the It- 
matrix calculation in this work. The black dots indicate the experimental 
asymmetry parameter of Southworth et al.— averaged over these five NO + 
electronic states 
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FIG. 1. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
N+X 2 £+ state. The cross sections are shown in unit of mega barn (Mb), equal to 10 22 m 2 . Our 
CASSCF target results and SCF target results are represented as thick full lines and thick dashed 
lines, respectively. The previous theoretical results of Montuoro and Moccia^ and Stratmann et 
al.~ are shown as thin dotted and thin dashed lines, respectively. Experimental data included in 
the panel (a) are taken from Hamnett et al.— , Samson et al.— and Plummer et al^. Experimental 
data in the panel (b) are taken from Marr et al4^ and Southworth et al4£ 
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FIG. 2. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
N^A 2 !! state. The other details are the same as in Fig. [TJ 
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FIG. 3. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
N^i? 2 S+ state. The other details are the same as in Fig. [TJ 
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FIG. 4. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
N^"2 2 S^ state. The symbol represents the cross section for the "Z" state of N^" in Hamnett et al.— 
The other details are the same as in Fig. [TJ 
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FIG. 5. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
NO+X 1 S+ state. Our CASSCF target results and SCF target results are represented as thick full 
lines and thick dashed lines, respectively. The previous theoretical results of Stratmann et al*2£ 
are shown as thin dashed lines. Experimental data in the figure are taken from Southworth et al4£ 
and Iida et al.— 
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FIG. 6. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
NO + 6 3 n state. The other details are the same as in Fig. [5j 
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FIG. 7. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
NO+A 1 !! state. The other details are the same as in Fig. [5j 
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FIG. 8. Photoionization cross section (a) and asymmetry parameter (b) for ionization of the 
NO + ti; 3 A, a 3 S + , A 1 H~ , 6' 3 S~ and P^ 1 A states, obtained by the R-matrix calculation in this 
work. The black dots indicate the experimental asymmetry parameter of Southworth et al4£ 
averaged over these five NO + electronic states. 
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TABLE I. Vertical ionization potentials of N2 molecule for the lower N^" ionic states (eV). Exper- 
imental values are taken from Baltzer et al£& 
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TABLE II. Ionization potentials of NO molecule for the lower NO + ionic states (eV). Our IPs are 
calculated for the vertical ionizations, while the experimental IPs taken from Albtitton et al.— are 
for the adiabatic ionizations. 
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